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Abstract 

The purpose of the present paper is to study the neutrino oscillations as they may 
appear in the low energy neutrinos emitted in triton decay: 

\H ->| He + e~ + v e 

with maximum neutrino energy of 18.6 KeV . Such low energy neutrino oscillations 
can only be detected via neutrino electron scattering. The novel feature is the fact 
that the cross section is sensitive to both charged and neutral currents. Thus one can 
simultaneously probe the electron neutrino disappearance as well as the appearance 
of the other two flavors. By folding the differential cross section with the neutrino 
spectrum one can study the dependence of the oscillation parameters as a function 
of the electron energy. The technical challenges to this end can be summarized 
as building a very large TPC capable of detecting low energy recoils, down to a 
few 100 eV, within the required low background constraints have been previously 
described [1]. We only mention here that the oscillation length involving the small 
angle S13 = sin #13, directly measured in this experiment, is fully contained inside 
the detector. A sensitivity of a few percent for the measurement of the above angle 
will be achieved. 
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1 Introduction. 



The discover of neutrino oscillations can be considered as one 
of the greatest triumphs of modern physics. It began with at- 
mospheric neutrino oscillations [2] interpreted as — > v T oscil- 
lations, as well as v e disappearance in solar neutrinos [3]. These 
results have been recently confirmed by the KamLAND exper- 
iment [4], which exhibits evidence for reactor antineutrino dis- 
appearance. As a result of these experiments we have a pretty 
good idea of the neutrino mixing matrix and the two indepen- 
dent quantities Am 2 , e.g ml — ml ans m 2 — m 2 |. Fortunately these 
two Am 2 are vastly different 

Am 2 (21)m 2 21 = \m\ - m\\ = (5.0 - 7.5) x 10 _5 (eF) 2 

and 

Amj 2 = \m\ - m\\ = 2.5 x Kn 3 (eT0 2 . 

This means that the associated oscillation lengths are very dif- 
ferent and the experimental results can effectively be described as 
two generation oscillations. This has been followed by an avalanche 
of interesting three generation analyses [5] -[9]. 

In all of these analyses the oscillation length is much larger than 
the size of the detector. So one is able to see the effect, if the 
detector is placed in the right distance from the source. It is, 
however, possible to design an experiment with an oscillation 
length of the order of the size of the detector. In this case one can 
start with zero oscillation near the source, proceed to maximum 
oscillation near the middle of the detector and end up again with 
no oscillation on the other end. This is achieved, if one considers a 
neutrino source with the lowest possible average neutrino energy 
such as a triton source with a maximum energy of 18.6 keV. Thus 
the average oscillation length is 6.5m, which is smaller than the 
radius of 10m of a spherical TPC detector (for a description of 
the apparatus see our earlier work [1]. 

Since the neutrino energies are so small, the only possible detec- 
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tor sensitive to oscillations is one, which is capable of detecting 
electrons. With this detector can detect: 

• electrons which are produced by electron neutrinos via the 
charged current interaction, which are depleted by electron 
neutrino disappearance 

• electrons which may arise from the other two neutrino fla- 
vors due to the neutral current interaction. These two flavors 
are produced by the muon and tau neutrino appearance, 
which operates simultaneously. 

Since these two mechanisms are competing with each other and 
the subsequent electron cross section depends on the electron 
energy, one has a novel feature, i.e. the dependence of the ef- 
fective oscillation probability on the electron energy. Thus the 
results may appear as disappearance oscillation in some kinemat- 
ical regime and as appearance oscillation in some other regime. 
In this paper we will examine how these features can best be 
exploited in determining the neutrino oscillation parameters. 



2 Elastic Neutrino Electron Scattering 



For low energy neutrinos the historic process neutrino-electron 
scattering [10] [11] is very useful. The differential cross section 
[12] takes the form 

da fda\ ^ fda\ (2 1) 

dT \dT/ wea k \dT J em 

The cross section due to weak interaction alone becomes: 

{%) =^[(9v + 9A? (2.2) 

\al / weak 27T 

+ (9v-9A?[l-^f + (91-91)^} 
gy = 2 sin 2 6 W + 1/2 [y e ) , g v = 2 sin 2 6 W - 1/2 (v^ v T ) 



3 



g A = 1/2 (v e ) , g A = -l/2 {y^v T ) 
For antineutrinos g A —> — gA- To set the scale 



GpTUe 

2tt 



= 4.45 x 10 



-48 



cm 



keV 



(2.3) 



The electron energy depends on the neutrino energy and the scat- 
tering angle and is given by: 



T = 



2 m e {E v cos 6)' 



{m e + E v f-{E v cos Of 
The last equation can be simplified as follows: 



T 



m f 



The maximum electron energy depends on the neutrino energy. 
For E v = 18.6 KeV one finds that the maximum electron kinetic 
energy is: 



T 



1.27 KeV 



I what follows, whenever appropriate, we are going to average 
our results with the neutrino spectrum shown in Fig. 2.1. 



0.08 



. 04 




7.5 10 12.5 15 17.5 



Fig. 2.1. The neutrino spectrum distribution function for the triton source (normal- 
ized to unity). 
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3 NEUTRINO OSCILLATIONS IN THE NOSTOS EXPERIMENT 



The electron neutrino is given [1] by: 

C Vi + S V 2 + Sl3 ^3 



V e = 



I+S13 



where s « sin6 so i ar c « cos6 so i ar) 9 so i ar determined from the solar 
neutrino data, and S13 = sinOis is a small quantity constrained 
by the CHOOZE experiment to be small. The electron disappear- 
ance oscillation probability is given by: 

P(z/ e -> v e ) = 



1-4 



s 2 c 2 sin 2 (7r^) + (s 2 sin 2 (7r^) + c 2 sin 2 (7r^)j s 2 3 



(1 + s\,f 
Assuming Am^ « Am^ we find: 



1 - 



P(> e -> i/ e ) « 
(sin 2^ ar ) 2 sin 2 (7r^) + 4s 2 3 sin 2 (7r^) 



(1 + sf 3 ) 2 

where the term proportional to s 2 3 is relevant for NOSTOS (L21 = 
5OL32). The NOSTOS experiment, in addition to the electronic 
neutrinos, which are depleted by the oscillations, is sensitive to 
the other two neutrino flavors, which can also scatter electrons 
via the neutral current interactions. These flavors are generated 
via the appearance oscillation: 

a#e 

(sin2^ ar ) 2 sin 2 (7r^) + 4s 2 3 sin 2 ^ 
(l + * 2 3 ) 2 

Thus the number of the scattered electrons, which bear this 
rather unusual oscillation pattern, is proportional to the (v e , e ) 
scattering cross section with a proportionality constant given by: 
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P(i/ e - i/ e ) « 1 T) 

(sin2^ ar ) 2 sin 2 (7r^) + 4s 2 3 sin 2 (7r^) 

(l + s? 3 ) 2 

(3.4) 

where x(E v ,T) is one minus the fraction of the cross section of 
flavor or v T divided by that of u e , i.e. 



X (E V , T) = 4sin 2 6 w [2E 2 V - ATE V + T(2T - m e ) / 
[£ 2 (1 + 4 sin 2 9 W + 8 sin 4 W ) - 2E V T{\ + 2 sin 2 ^ 
+ (1 + 2 sin 2 9 W )T(2(T - m e ) sin 2 6 W + T)] 



(3.5) 



The function x{E v , T) in the region: 



Tfl ! / 777 ^ / 

E v + -±- -y/8m e E v + m 2 < T < + + -y/8m e E„ 



+ 777/ 



becomes negative. This means that the cross section due to 
or v T becomes larger than that of v e . This can also be seen by 
plotting the relevant cross sections, see Figs 3.2 and 3.3. We see 
that the cross section due to flavor or v T is rising, while the 
one due to v e is dropping, see Fig. 3.3. This means that the 




Fig. 3.2. The differential cross (v e ,e~) section, in arbitrary units, as a function of 
the neutrino energy for electron energies from left to right 0.2, 0.6 and 1.0 keV. The 
the dashed curve corresponds to the cross section due to u a , a = \i or r. 

function x(Ey, T) eventually changes sign, see Fig. 3.4. If this oc- 
curs in a kinematically interesting region, it will lead to neutrino 
oscillation results different than those naively expected. 
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Fig. 3.3. The differential cross (z/ e ,e _ ) section , in arbitrary units, as a function of 
the electron energy for neutrino energies from left to right 7.6, 13.1 and 18.6 keV. 
The the dashed curve corresponds to to the cross section due to u a , a = \i or r. 




Fig. 3.4. On left the function x(Ev, T) is shown as a function of the electron energy 
for neutrino energies 7.6 keV (thin line), 13.1 keV (dashed line) and 18.6 KeV (thick 
line) . On the right we show the same quantity as a of the neutrino energy for electron 
energies 0.2, 0.4, 0.6, 0.8, 1.0, 1.2 keV increasing to the right. 

4 VARIOUS PATTERNS OF OSCILLATION PROBABILITIES 

The oscillation probabilities for NOSTOS, i.e. near the source, as 
well as experiments further from the source for sin 2 29i3 = 0.17 
with and without the function x{E v , T) are shown in Figs 4.5-4.8. 
For sin 2 26 13 = 0.085 in Figs 4.9-4.12 and for sin 2 26 13 = 0.045 in 
Figs 4.13-4.16 




Fig. 4.5. The oscillation probability at distances larger than the NOSTOS size with- 
out the function x{E v ,T). One can see the oscillation associated with Amf 2 - If 
sensitive to to sm 2 2# 13 = 0.17 one will also see on top of this the small oscillation 
associated with Am|. From left to right, E v = 7.6, 13.1 and 18.6 keV. 
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Fig. 4.6. The same as in Fig. 4.5 after including the function x(E u ,T). The results 
now depend on the electron energy. For E v = 7.6 (left) we show the results for 
T = 0.200, 0.213 and 0.220 keV indicated by thin, dashed and thick line respectively. 
For E v = 13.1 (middle), T = 0.200, 0.419 and 0.639 keV, while for E„ = 18.6 (right), 
T = 0.200,0.731 and 1.262 keV. 




Fig. 4.7. The oscillation probability expected in NOSTOS with sin 2 29i3 = 0.17 
corresponding to, from left to right, E u = 7.6,13.1 and 18.6 keV, without the 
function x(E u , T). 




Fig. 4.8. The same as in Fig. 4.7 after including the function x{E v ,T) and the 
neutrino and electron energies as indicated in Fig 4.6. 




Fig. 4.9. The same as in Fig. 4.5 with sin 2 2d 13 = 0.085. 




Fig. 4.10. The same as in Fig. 4.6 with sin 2 29 13 = 0.085. 
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Fig. 4.11. The same as in Fig. 4.7 with sin 2 29 13 = 0.085. 




Fig. 4.12. The same as in Fig. 4.8 with sin 2 29 13 = 0.085. 




Fig. 4.13. The same as in Fig. 4.5 with sin 2 29 13 = 0.045. 




Fig. 4.14. The same as in Fig. 4.6 with sin 2 29 13 = 0.045. 




Fig. 4.15. The same as in Fig. 4.7 with sin 2 29i 3 = 0.045. 




Fig. 4.16. The same as in Fig. 4.8 with sin 2 29i 3 = 0.045. 
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5 AVERAGED NEUTRINO OSCILLATION PROBABILITIES 



With the above ingredients one can obtain an effective neutrino 
oscillation probability for various electron energies by averaging 
the expression of the oscillation probability, Eq. (3.4), over the 
weak differential cross section, Eq. 2.3, and the neutrino spec- 
trum, see Fig. 2.1, i.e. 




Fig. 5.17. The average oscillation probability as a function of the radial distance 
L (in meters). We have integrated over the neutrino spectrum, but the function 
X{E U ,T) was not included. From left to right sin 2 26 n = 0.045,0.085,0.170. The 
depth is increasing with T 

Oscillations for sin 2 26n = 0.17 as seen in distances far greater 
than the NOSTOS experiment as well as in the NOSTOS ex- 
periment are shown in Figs 4.5-4.8 both with and without the 
function x(E v ,T). Results for sin 2 26ys = 0.085,0.045 are shown 
in Figs 4.9-4.12 and 4.13-4.16 respectively. 




Fig. 5.18. The average oscillation probability as a function of the radial distance L. 
We have integrated over the neutrino spectrum and included the function x{E v , T). 
From left to right sin 2 26i 3 = 0.045,0.085,0.170. Labeling as in Fig. 5.17. 
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6 RESULTS 



We now are going to present what the NOSTOS experiment is 
going to measure, namely the differential cross section averaged 
over the neutrino spectrum, i.e. 

= ^p(u e ^u e )y^f e (T)y (6.7) 

Or 

-< % = * P("e "> Ve) ^< fe(T) >- (6.8) 

-< fe(T) y= J {%) MEjdE,, (6.9) 

J \al J weak 

The results obtained are shown in Figs 6.19-6.20. We see that, 
unfortunately, the unexpected oscillation features are washed out. 
They stood out mainly due to the fact that we normalized the 
amplitude in the usual way, i.e. so that the oscillation length in- 
dependent part of the amplitude is equal to unity. They happen, 
however, to occur in a kinematically unfavored region, where the 
cross sections is small. Indeed the differential cross section, av- 
eraged over the neutrino spectrum, as we have seen can also be 
obtained by multiplying the average oscillation probability with 
the function given by Eq. 6.9 (see Eqs 6.8 and 5.6). The latter 
quantity is shown in Fig. 6.22. We thus see that the contribu- 
tion of the high energy electrons becomes very small. One can, 
of course, obtain the total cross section by integrating the differ- 
ential cross section over the electron energy as a function of the 
distance from the source. The thus obtained results are shown in 
Fig. 6.23. 

Experimentally, of course, one measures events. The number of 
events is obtained by: 
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Fig. 6.19. The differential cross (u e ,e ) section,^ ^> m units of 

= 4.5 x 10~ 52 ^7, as a function of the source-detector distance averaged 
over the neutrino energy for electron energies from top to bottom and left to right 
0.2, 0.4, 0.6, 0.8, 1.0 and 1.2 keV. The results shown correspond to sin 2 29iz = 0.170 

• Inserting in the relevant cross sections the scale of week inter- 
action 



= 4.50 x 10' 



2tt keV 

• Considering the total number of neutrinos emitted by the source 
of 20 Kg of tritium 

N v = — = 4.00 x 10 27 

3 x 1.67 x 10- 27 

for any given time t the number of neutrinos emitted must 
be multiplied by the fraction 1 — e~ t ' T with r = ~ 17.7?/ 
(^1/2 = 12.3?/). This yields the neutrino flux: 

• Assuming that in the target the number of electrons per unit 
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Fig. 6.20. The same as in Fig. 6.19 for sin 2 29 13 = 0.085 

volume is 

ry P a a m27 -3 P Z 300 

n e = Z- — = 4.4 x lCTra 6 — 

kT 10 Atm 18 T 

[Z the atomic number, P is the pressure and Tq the tempera- 
ture). 

• The total number of events with energies between T and T+dT 
in the length region between L and L + dL is 

diV = ^ ^ n e 4ttL 2 dTdL = N v ^= n e dTdL. 
dT dT 

For any given time t we get only the fraction 1 — e _ ^ T of these 
events, e.g 5% the first year, 11% for two years of running, 16% 
after 3 years and 43% after 10 years. 

The thus obtained results for the differential rate d ^ N dL , during 
the life time of the source (20 Kg of tritium), employing |§Ar at 
P = 10 Atm, are shown in Fig. 6.24-6.26. The total rate per 
unit length, obtained by integrating the differential rate of Fig. 
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Fig. 6.21. The same as in Fig. 6.19 for sin 2 26 13 = 0.045 



f 




0.2 0.4 0.6 0.8 1.2 

Fig. 6.22. (/ e (T)), the differential cross (u e , e~) section without oscillations, in units 

of ° F ^ £ = 4.5 x 10~ 52 ^j, averaged over the neutrino spectrum, as a function of 
the outgoing electron energy in keV. 



6.24-6.26 with respect to the electron energy, is shown in Fig. 
6.27. 
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Fig. 6.23. The total cross section, in units of = 4 5 x \g 52 C m 2 , as a func- 

tion of the source detector distance. From top to bottom and left to right for 
sin 2 29 13 = 0.170, 0.85, 0.045. 

7 Discussion and Conclusions 

In the experiment one will measure the differential rate, i.e. the 
number of electrons detected per unit energy and per unit length, 
as a function of the distance L from the source, i.e. the radius 
of the sphere. Integrating the differential rate one obtains the 
number of electrons per unit length as a function of L. From Figs 
6.24-6.26 and 6.27, obtained with a source of 20 Kg of tritium over 
its lifetime and using Ar at 10 Atm, we see that the L-dependence 
of the event rates resembles closely that of the corresponding 
cross sections shown in Figs 6.19-6.21 and 6.23. To get the event 
rate after some time t of running, we must multiply by the factor 
1 — e~*/ T , t = 17. 7y, e.g. 5% during the first year etc. The event 
rates scale with the atomic number Z . From these plots we see 
that: 

• In the case of the differential rate: 

(1) At high electron energies it increases as a function of L, i.e. it 
behaves like appearance oscillation. The electrons are mainly 
produced by the other two appearing neutrino flavors. This 
is a novel feature of the neutrino electron scattering. 

(2) At low electron energies, where the cross section is large, one 
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Fig. 6.24. The differential rate d ^ dL (per keV-meter) for Ar at 10 Atm with 20 Kg 
of tritium as a function of the source-detector distance (in m ), averaged over the 
neutrino energy, for electron energies from top to bottom and left to right 0.2, 0.4, 
0.6, 0.8, 1.0 and 1.2 keV. The results shown correspond to sin 2 29 13 = 0.170. This 
rate must be multiplied by 1 — e~ l l T to get the number of events after running time 
t. Thus only 5% of these are expected during the first year of running. 

recovers the usual features of electron flavor disappearance. 

(3) The difference between the maximum and the minimum is 
at the level of 8,4 and 2 % for sin 2 26 13 = 0.170, 0.085 and 
0.045 respectively. 

(4) The location of the maximum depends on the electron energy. 
• In the case of the total cross section: 

(1) The percentage of the difference between the maximum and 
the minimum is about the same as before. 

(2) The maximum of the oscillation probability occurs at L = 7.5 
m. 

It is clear that with such a large number of events the above fea- 
tures of neutrino oscillation can be exploited to measure sm 2 2/9i3 
and perhaps to accurately determine the oscillation length. 
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Fig. 6.25. The same as in Fig. 6.24 for sin 2 26 13 = 0.085. 
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